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ABSTRACT 

Detailed studies of the Spectral Energy Distribution (SED) of normal galaxies are been increasingly used in order to understand the 
physical mechanism dominating their integrated emission, mainly due to the availability of high quality multiwavelength data from the 
UV to the far-infrared (FIR). However, systems hosting dust-enshrouded nuclear starbursts and/or an active galactic nucleus (AGN) 
due to an accreting supermassive black hole, are especially challenging to study. This is due to the complex interplay between the 
heating by massive stars and the AGN, the absorption and emission of radiation from dust, as well as the presence of the underlying 
old stellar population. 

We use the latest release of CIGALE, a fast state-of-the-art galaxy SED fitting model relying on energy balance, to study in a self 
consistent manner the influence of an AGN in estimating both the star formation rate (SFR) and stellar mass in galaxies, as well as 
to calculate the contribution of the AGN to the power output of the host. Using the semi analytical galaxy formation model galform, 
we create a suite of mock galaxy SEDs using realistic star formation histories (SFH). We also add an AGN of Type 1, Type 2, or 
intermediate type whose contribution to the bolometric luminosity can be variable. 

We perform an SED fitting of these catalogues with CIGALE assuming three different SEHs: a single-exponentially-decreasing 
(Ir-dec), a double-exponentially-decreasing (2T-dec), and a delayed SEH. 

Constraining the overall contribution of an AGN to the total infrared luminosity {frac/^Q^) is very challenging for frac/^Q^ <20%, 
with uncertainties of ~5-30% for higher fractions depending on the AGN type, while FIR and sub-mm are essential. The AGN power 
has an impact on the estimation of M, in Type 1 and intermediate type AGNs but has no effect for galaxies hosting Type 2 AGNs. We 
find that in the absence of AGN emission, the best estimates of M, are obtained using the 2T-dec model but at the expense of realistic 
ages of the stellar population. The delayed SFH model provides good estimates of M, and SFR, with a maximum offset of 10% as 
well as better estimates of the age. 

Our analysis shows that the underestimation of the SFR increases with frac/^Q^ for Type 1 systems, as well as for low contributions 
of an intermediate AGN type, but it is quite insensitive to the emission of Type 2 AGNs up to fracACN ~45%. 

A lack of sampling the FIR, or submm domain yields to a systematic overestimation of the SFR (<20%), independent from the 
contribution of the AGN. Similarly the UV emission is critical in accurately retrieving the M, for Type 1 and intermediate type AGN, 
and the SFR of all of the three AGN types. 

We show that the presence of AGN emission introduces a scatter to the SFR-M, main sequence relation derived from SED fitting, 
which is driven by the uncertainties on M,. 

Einally, we use our mock catalogues to test the popular IR SED fitting code DecompIR and show that fracACN is underestimated but 
that the SFR is well recovered for Type 1 and intermediate types of AGN. The /rac^Civ, SFR, and Ljf, estimates of Type 2 AGNs are 
more problematic due to a FIR emission disagreement between predicted and observed models. 

Key words. Galaxies: fundamental parameters, active 


1. Introduction 

The formation of galaxies and their evolution with cosmic 
time are open problems in current astrophysical research. 
Understanding the assembly of galaxies and the built-up of their 
stellar populations is challenging because the relevant phys¬ 
ical processes are complex, interconnected and operate on a 


large range of scales. Gas inflows from the intergalactic medium 
(IGM) for example, are important for supplying galaxies with 
fresh material that can be turned into stars or feed supermassive 
black holes (SMBH) at their centers. Feedback processes asso¬ 
ciated with stellar evolution of Active Galactic Nuclei (AGN) 
can drastically modify the physical conditions of the Inter-Stellar 
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Medium (ISM) thereby, affecting the formation of new stars. 
Supernovae explosions enrich with heavy metals the gaseous 
component of galaxies and their environments and therefore 
change the composition of the ISM with implications on star- 
formation. The density of IGM on large scales or interactions 
with nearby galaxies also have an impact on the ISM of individ¬ 
ual systems and can modify their evolutionary path. 

One approach for shedding light into the physics of galaxy 
formation and evolution is population studies. Multi-wavelength 
observations provide information on galaxy properties such as 
stellar mass, star-formation history (SFH), gas content, AGN ac¬ 
tivity, kinematics, structural parameters or position on the cos¬ 
mic web. Each of these observationally determined parameters 
probe different physical processes. Exploring correlations be¬ 
tween them for large samples can therefore provides insights 
on how galaxies form and evolve in the Universe. Eor exam¬ 
ple, the discovery of the bimodality of galaxies on the colour 
(proxy to star-formation history) versus stellar-mass plane (e.g. 
Strateva et al.T20ml [Baidry et al.|2004HBell et al.|2004l|Baldry| 


et al.|2006|l has b^n interpreted as evidence for star-formation 


quenching that may be driven by either internal (e.g. feedback) 
or external (e.g. environment) processes (e.g. [Bell et al.||2004| 
|Eaber et al.||2()07t [Schawinski et al.||2014[ ). The star-formation 
main sequence (MS), i.e. the relatively tight correlation between 
stellar mass and star-formation rate (SER) of star-forming galax- 


The overall consensus is that stellar masses can be reliably 
constrained, although systematics up to the 0.5 dex level remain, 
depending on the adopted Initial Mass Eunction, stellar popu¬ 
lation libraries, functional form of the model SEH and the im¬ 
plementation of dust attenuation. The SER of galaxies is often 
determined using as tracers the observed luminosity at specific 
wavelengths (e.g. |Kennicutt|19^|Hopkins|2004| l. Nevertheless, 
studies that fit the full SED to derive SER find an overall rea¬ 
sonable agreement with estimates based on specific tracers (e.g. 
ISalim et al.|2007||200^[Walcher et al.|20081 l. 

One aspect of galaxy evolution that has developed consid¬ 
erably in recent years is the relation to SMBH growth. This 
has been motivated by the tight correlations between proxies 
of the stellar mass of local spheroids and the mass of the black 
hole at their centres (e.g. Kormendy & Ho|2013 and references 
therein). One interpretation of these correlations is that of the 
co-evolution of AGN and galaxies dictated by a common gas 
reservoir that forms new stars and also feeds the central black 
hole. An alternative explanation is that of AGN outflows that 
heat and/or expel the cold gas component of galaxies thereby, 
regulating the formation of new stars and ultimately the accre¬ 
tion onto the central SMBH itself (e.g[Silk & Rees|1998[|Eabian| 


1999t King|20d3 2005 [ Di Matteo et al.|2005| l. There is indeed 


increasing observational evidence for AGN-driven winds at both 
the local Universe and high redshift (e.g 


ies (e.g. [Salim et al. 

2007} 

Noeske et al. 

2007b| Elbaz et al. 

Blustin et al.|2005 

Tombesi et al.|2010 

2012 [Saez & Chartas 

2007 [ Rodighiero et a 

T20T1 

Speagle et al. 

201411, suggests that 

2011 

[Lanzuisi et a 

.|[2012l [Cicone et a 

20121 2014 Harrison 


the bulk of the stars in the Universe form via secular processes 
rather than in violent events, such as mergers. Addtionally, the 
slope and redshift evolution of the MS normalization has been 
discussed in context of feedback processes and gas exhaustion 
with cosmic time, respectively (e.g. Noeske et al.|2007a{|Zheng 
|et al.||2007] [Tacconi et al.||20l3] |Guo et akpoisf . Corrdations 

between galaxy structural parameters, SER and stellar mass (e.g. 
Kauffmann et al.|2004t|Wuyts et al. 201 1[ Bell et al.|2012[ Lang 


et al.|2014| l suggest that the formation of galaxy bulges and the 

quenching of star formation are likely related to the same under¬ 
lying processes. 

Among the different galaxy properties that are accessible to 
observations, the stellar mass and SER play an important role 
in galaxy evolution studies. This is not surprising. Both quan¬ 
tities provide a measure of the SEH of the galaxy, either inte¬ 
grated over its lifetime (stellar mass) or averaged in the last few 
tens to few hundred million years (instantaneous SER). Indeed, 
all galaxy properties show strong correlations with either stel¬ 
lar mass, SER or both (e.g., [Kauffmann et al.||2004| l. This has 
led to an increasing refinement in methodology to provide more 
accurate and less biased observational constraints to the stellar 
mass and SER of individual galaxies. Among the different ap¬ 
proaches, the one that has been extensively used in the literature 


mon these outflows are and if they are energetically important to 
affect the ISM of their hosts. One way to approach these issues 
and place AGNs in the context of galaxy evolution is statisti¬ 
cal studies of the host galaxy properties of large AGN samples. 
Questions that this approach could address include, when during 
the lifetime of a galaxy accretion onto the central SMBH is trig¬ 
gered, how black hole growth is related to star-formation, and 
whether AGNs affect their host galaxy properties. 

Multi-wavelength survey programs in the last decade have 
started addressing these questions. Ear-infrared (EIR) and sub- 
mm surveys with Herschel for example, measure the mean SER 
of AGN hosts via stacking methods and show that AGN, on aver¬ 
age, lie on the main sequence star-formation at all redshifts (e.g. 
Mullaney et al.||2()12t jSantini et al. 2012 Rosario et al.||201^ 


2013|). This finding suggests that at least in an average sense the 


is the use of stellar population synthesis models (e.g. Bruzual 
|& Charlot||2003| |Marakon]|2005| l to generate Spectral Energy 
Distribution (SED) templates for different star-formation histo¬ 
ries and then fit them to the broad-band photometry of galax¬ 
ies (e.g. jWalcher et al.j 2011| and references therein). This is 
driven by the explosion in recent years in the quality and quan¬ 
tity of multi-wavelength imaging surveys, which provide UV to 
infrared (IR) SEDs of large galaxy samples over a wide range 
of redshifts. This approach for inferring galaxy parameters has 
also been extensively tested to identify limitations and potential 
sources of systematics (e.g. Pozzetti et al. [20071 Marchesini et al.j 


same physical processes govern the formation of stars in galax 
ies and the growth of black holes at their centers (Silverman 
[et al.[|2009] [Georgakakis et aL][2011[ [Mullaney et al. ]2012[ l. 
Studies of the star-formation properties of individual AGN hosts 
(rather than averages over populations) suggest a wide range 
of SEH. The rest-frame colors of AGN hosts ( [Aird et al.|2012[ 
Georgakakis et al.|2014 Azadi et al.|2014|l or their position on 


the color-magnitude diagram (e.g.|S4nchez et al.|2004[ Nandra 
[et al.|2007[ [Schawinski et ar][2009[ ) indicate that they are scat¬ 
tered at the red sequence of passive galaxies, the star-forming 
cloud and the green valley in between. A higher incidence of 
AGN among green valley galaxies is claimed ( [Nandra et al.| 
2007 Schawinski et al. 2009[), i.e. systems with colors in the 


20091|Conroy et al.|2009 

Wuyts et al.|2009( Ilbert et al. 

2010 

Michalowski et al.[[2012[ 

Pforr et al.[[2012[ Banerji et al. 

2013 


Schaerer et al.|2013[[Mitchell et al.|2013[[Buat et al.|2014[ ). 


transition region between the red sequence and the blue cloud. 
This can be interpreted as evidence for AGN being responsible 
for the quenching of star-formation in galaxies and hence, their 
transition from blue star-forming to red and dead systems. At the 
same time however, the importance of constructing appropriate 
control samples when determining the fraction of AGN among 
galaxies is also emphasized. The evidence for an increased AGN 
fraction among green valley galaxies is less strong, or may even 
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disappear once the stellar mass of AGN hosts and the control 


galaxy samples are matched (e.g. Silverman et al. 2009 Xue 
et al.|201^ Aird et ar||2012| l. Constraints on the ste lar ages of 
AGN hosts further suggests that black hole accretion preferen¬ 
tially occurs few hundred to few thousand million years after the 
peak of star-formation (Kauffmann et al. 2003 a[ Wild et al.|2010[ 


Hernan-Caballero et al. 20141, i.e. during a period when galaxies 
can be identified as post-starbursts (e.g. Goto|200^ Georgakakis 


et al.|200^. A time lag between the peak of star-formation and 


the black hole growth may pose a problem to AGN-driven star- 
formation quenching scenarios. 

The studies above demonstrate the importance of measur¬ 
ing the star-formation and stellar mass for AGN hosts. This is 
not only for understanding the relation between star-formation 
history and black hole growth in individual systems, but also 
for constructing appropriate control samples of inactive (non- 
AGN) galaxies. The determination of stellar masses or SFR for 
AGN hosts is challenging. Emission from the central engine con¬ 
taminates or may even dominate the underlying galaxy light, 
thereby rendering the determination of galaxy properties, via e.g. 
SED fitting methods, difficult. Different approaches have been 
adopted in the literature to address this issue. Standard SED fit¬ 
ting methods are often applied to determine host galaxy proper¬ 
ties only for low luminosity and obscured AGN, under the rea¬ 
sonable assumption that contamination in these objects is small. 
A potential problem with this approach is that under certain 
models for the co-evolution of AGN and galaxies (e.g |Hopkins| 
|et al.|[2006l l different levels of obscuration and different accre¬ 
tion luminosities correspond to different evolutionary phases of 
the black hole growth. The selection against unobscured and lu¬ 
minous AGN might introduce biases and lead to erroneous con¬ 
clusions on the relation between black hole accretion and galaxy 
formation. An alternative approach is to add AGN templates to 
the stellar template library and perform SED deconvolution to 
separate the galaxy and AGN components (e.g. Bongiorno et al. 
|2012[ [Lusso et al.||201l| |2013[ [Rovilos et al.| 2014[ ). This can 


be powerful, although degenerations between AGN and galaxy 
templates may introduce systematics in the determination of host 
galaxy parameters. 

To date there are still few studies that explore and quantify 
the impact of AGN contamination on the determination of the 
underlying host galaxy properties via SED fitting (|Wuyts et al. 


2011 1 Hayward & Smith||2014 1 . In this paper, we address this 


using the bayesian-based SED fitting code CIGALE (|Noll et al? 


2009 Burgarella et al., in prep; Boquien et al., in prep). To do 
so, we follow the method developed in Mitchell et al.|(2OT3]l and 


use simulated galaxies from the SAM (Semi-Analytica Model) 
code GALFORM from which we know the exact value of the stellar 
mass and SFR. After building the UV-to-submm SEDs of the 
GALFORM objects, and adding an AGN contribution, we will study 
our ability to retrieve the original properties using CIGALE. 

The paper is organized as follows. Eirst, we introduce gal- 
form and present the derived star formation history (SEH) of a 
hundred galaxies at z = 1 (Section |^. We describe how we 
build our mocks SEDs using the modeling function of CIGALE 
in Section]^ and how we perform the SED fitting in Section]^ 
In the same section, we also present the comparison between 
the true values of the stellar masses and SFRs and the outputs 
of CIGALE in normal galaxies as well as in AGN host galax¬ 
ies. Einally, we discuss the impact of our results on the SER-M, 
relation in Section In Section we use our mock samples 
to evaluate the performance of the popular IR SED fitting code 
DecompIR ( |Mullaney et al. [201 l| l. The purpose of this work is 
not to test thoroughly the validity of AGN emission models, nei¬ 


ther the ability of SED fitting codes to accurately retrieve the 
properties of the AGN through its IR emission but to analyze the 
impact of this emission on the derivation of the basic host galaxy 
properties. 

In this work, we assume that Q.„, - 0.25, Qa = 0.75, and 
Ho - 73 km s“^ Mpc^'. These values are used because the gal- 
form model is built on top of merger trees generated from the 
original Millenium simulation ( [Springel et al.|2005 i, which as¬ 
sumes a WMAPl like cosmology. This choice does not impact 
our results, as we only use this cosmology to convert physical 
quantities to observables, without comparisons with other works 
based on more recent cosmological parameters. All of the stellar 
masses and SFRs are provided assuming a IMF of Salpeter. 


2. Building a realistic mock galaxy sample 

In this paper we explore the problem of AGN/stellar light de¬ 
composition via SED template fits to multi-wavelength broad¬ 
band photometric data. The goal is to quantify how well the 
galaxy properties of AGN hosts, such as stellar mass and star- 
formation rate, can be measured from broad-band photometry. 
Our approach is to use mock galaxies extracted from the gal- 
form SAM ( |Cole et al.|2000| l, for which the SFHs, stellar masses 
(Mt) and instantaneous SFR are known. This information allows 
the construction of UV to submm SEDs for the mock galaxies. 
AGN templates are then added to the galaxy emission to gen¬ 
erate composite SEDs. These are then integrated within broad¬ 
band filters to generated mock photometric catalogues. The fit¬ 
ting modules of the CIGALE code are then applied to the mock 
galaxy photometry to decompose AGN from stellar light. Each 
of these steps are described in the following sections. 


2.1. Simulations of realistic star formation histories 


To generate our grid of mock SFHs, we use the galform semi- 
analytic galaxy formation model to simulate the assembly of the 
galaxy population within the context of the ACDM model of 
structure formation ( |Cole et al.|2000[[Bower et al.|2006[[Benson] 
& Bower]|2010|l. The model is constructed on top of dark mat¬ 


ter halo merger trees extracted from the Millennium dark matter 
N-body simulation ( [Springel et al.|2005] l. Within each halo, the 
baryonic mass is divided into hot and cold gas along with stellar 
disk and bulge components. The model then solves a set of cou¬ 
pled differential equations that describe how mass is exchanged 
between these different components. Star formation in the model 
is split into quiescent star formation that occurs in galaxy disks 
and bursts of star formation that are triggered by galaxy merg¬ 
ers and disk instabilities ( Baugh et al.|2005] Bower et al.|2006 
[Lagos et al.|2011| . 

Star formation histories are extracted from the fiducial ver¬ 


sion of the model presented by Mitchell et al. (20141. We ran¬ 


domly select a total of 100 galaxies from that simulation at 
redshift z = 1. Motivated by observational evidence that AGN 
hosts are typically massive and lie, at least in average sense, on 
the main sequence of star-formation (e.g.,[Mullaney et al.|2012[ 


[Santini et al. [[2012] [Rosario et al.[[2012 [201 3| l, we also choose 

the mock galaxies to sit on the main sequence of the SAM at 
redshift z = 1 and to have high stellar masses. In particular they 
fulfill the following criteria: 


1. Specific Star-formation rates sSFR> 0.1 Gyr '. This cut sep¬ 
arates main star-formation sequence from passive galaxies in 
the model at z = 1 ( [Mitchell et al.|20l4| l. 

2. Stellar mass in the range 10 < log(M»/Mo) <11. 
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Fig. 2; Distribution of M, and SFR associated with the galform 
SFHs. The black lines show the mean value of the distributions. 


For these star formation histories, we sum over the star for¬ 
mation in all progenitors of the final galaxy and we also com¬ 
bine the stellar mass assembly of the disk and bulge together. 
Bursts of star formation in the model can occur over relatively 
short timescales in some cases and so we construct star forma¬ 
tion histories from the model to have high temporal resolution. 
We show ten examples of the produced SFFIs in Figure [T] The 
distribution of stellar mass and SFR at redshft z = 1 of the 
100 simulated galaxies selected for our analysis are presented in 
Figure]^ The mean values of the sample are logM,=10.37Mo 
and SFR=7.05M*.yr-'. 


2.2. Simulations of UV to submm SEDs 

CIGALI0 is a package that has two different and independent 
functions: an SED modeling function (Boquien et al., in prep) 
and an SED fitting function (Burgarella et al., in prep). The base¬ 
line functions of CIGALE are presented in [Noll et al. ( 2009) 1. 
Based on the same general principles as the original version 
of CIGALE, this entirely new version has been designed for a 
broader set of scientific applications and better performance. The 
latest version of the modeling function of CIGALE is briefly de¬ 
scribed in this section. 

The SED modeling function of CIGALE allows the con¬ 
struction of galaxy SEDs from the UV to the submm by assum¬ 
ing a stellar population library and star-formation histories pro¬ 
vided by the user. CIGALE builds the SED taking into account 
energy balance, i.e., the energy absorbed by dust in UV-optical 
is reemitted in the IR. 

In CIGALE, the SEH can be handled in two different ways. 
The first, is to model it using simple analytic functions (e.g. ex¬ 
ponential forms, delayed SEHs, etc). The second is to provide 
more complex (non-analytic) SEHs (e.g., |Boquien et al.|[20T4) l, 
such as those provided by the galform SAM. The stellar pop¬ 
ulation models of either [Maraston) ( |20051 l or |Bruzual & Chariot) 
( )2003) l are convolved with the adopted SEH to produce stellar 
SEDs, which are then attenuated by dust. The energy absorbed 
by the dust is reemitted in the IR using a choice of different dust 
templates ()Dale & Helou])2002) )Dale et al.))2014) )Draine & Li) 
[200^)Casey|2'0T^ 

CIGALE also allows the emission from AGN to be added to 




Eig. 3: SEDs generated by GIGALE using a particular SEH ex¬ 
tracted from GALFORM and adding on that different AGN tem¬ 
plates with different normalizations relative to the galaxy light. 
Each panel corresponds to one of the three Eritz et al. (2006) 
templates presented in Table the upper panel is for Type 1 
AGN, the middle panel is for the Intermediate AGN type, and the 
lower panel is for Type 2 AGN. SEDs are color-coded according 
to the contribution of the AGN to the total IR luminosity. The 
black solid lines at the bottom of each panel are the broad-band 
Alters of Table within which the model SEDs are integrated 
to generated mock photometric catalogues. The inset plots show 
the contribution of the three AGN components to the total SED 
for fracAGN—^0%. Red is for the dust torus emission, orange is 
the scattering component, and cyan the direct emission from the 
central AGN. 


et al. (2006 1 are adopted. These SEDs consist of two com¬ 


ponents. The first one is the isotropic emission of the central 
source, which is assumed to be point-like. This emission is a 
composition of power laws with variable indices in the wave¬ 


length range of 0.001-20/tm. The second component of the Fritz 
)et al.) ( )2006) l models is radiation from dust with a toroidal ge¬ 
ometry in the vicinity of the central engine. Part of the direct 
emission of the AGN is either absorbed by the toroidal obscurer 
and re-emitted at longer wavelength (1-1000 fim) or scattered by 
the same medium. Dust can be optically thick to its own radi¬ 
ation, thus requiring the numerical resolution of the radiative 
transfer problem. In )Fritz et SI) ( )2006) ) models, the conserva¬ 
tion of energy is always verified within 1% for typical solutions, 
and up to 10% in the case of very high optical depth and non¬ 
constant dust density. The choice of adding the )Fritz et al.) ( )200^ 
library into CIGALE is driven by the energy balance handling of 
the two components, which also matches the energy conserva¬ 
tion philosophy of CIGALE. Eurthermore, this library has been 


4 











































































Ciesla et al.: Derivation of stellar mass and SFR in AGN host galaxies 



Fig. 1: Examples of 10 SFHs of simulated galaxies extracted from the galform SAM. 


tested in numerous studies of the literature (e.g., Fritz et al. 2006 


Hatziminaoglou et al.|2008 2010t[Feltre et al.|2012|l. 

The relative normalization of these AGN components to the 
host galaxy SED is handled through a parameter which is the 
fraction of the total IR luminosity due to the AGN so that 


jAGN 

^IR 


■■ fracACN X 


TOT 


( 1 ) 


where is the AGN IR luminosity, fracACN is the con¬ 
tribution of the AGN to the total IR luminosity i.e. 

j^starburst _ Thus, estimating depends on the con¬ 

straints on fracACN- 

For each SFH provided by galform, we compute a host 
galaxy SED using th e |Maraston| (|2005|l ste llar population mod¬ 
els and assuming the Calzetti et al. (20001 extinction law. The 


amount of reddening, E{B - V)*, is chosen to be 0.2, motived 
by observational studies of AGN host galaxies ( |Kauffmann et al.| 
2003b Hainline et al.||2012| . The energy absorbed in UV-Opt- 


NIR is re-injected in IR, providing the normalization of the 
|Draine & Fi| ( |200'7] l models in order to maintain the energy bal¬ 
ance of the SED. The template is chosen following the results 
of|Magdis et al.| ( |2012[ ) who fit high-z galaxies with the Draine 
|& Li|P007^ models. The parameters used to model the galaxies 
are presented in Table [T] 


We select three AGN model templates from the Fritz et al. 
P006| l library to be added to the galaxy SEDs. These include a 
Type 1 AGN (i.e. unobscured), a Type 2 AGN (i.e. obscured) and 
a template that lies in-between the first two and is referred to as 
intermediate type. The latter model displays a power-law spec¬ 
tral shape in the mid-IR without strong UV/optical emission. 
Previous studies indicate that such an intermediate template is 
indeed necessary to represent the diversity of the observed SEDs 
of AGN ([Hatziminaoglou et al.|2008| 2009\ Feltre et al.|2012| l. 


The Fritz et al. ( 2006|l model parameters for the three templates 


are presented in Table 1 ■ The parameters for the Type 1 and 


Type 2 AGN templates are representative of local unobscured 
and obscured AGN, respectively ([Fritz et ar]|2006|l. For the in¬ 


termediate AGN type we adopt one of the [Fritz et al.j (|2006|l 
models with low equatorial optical depth ( Hatziminaoglou et al.[ 
2008[ 2009 [ Feltre et al. 2012[ Buat et al., in prep). The three 


AGN templates adopted in this paper are the minimum required 
to represent the variety of the broad-band AGN SEDs. There 
are strong degeneracies among different parameters of the Fritz 
[et al.[ ( [2006[ l library that cannot be broken by multiwavelength 
photometric data alone. Additional parameters combinations to 
those shown in Table [T] do not necessarily result in AGN broad¬ 
band SEDs that are distinctively different from our Type 1, Type 


^ We note an error in the definition of the angle relative to the line 
of sight if/ in [Fritz et al.H20(36t : i^=0° corresponds to a Type 2 AGN 
whereas an angle i^=90° is for a Type 1 (J. Fritz, private communica¬ 
tion). 


2 and intermediate type models. In Figure]^ we show the SEDs 
corresponding to one of the SFHs presented in Figure [T] 

The modeled SEDs are integrated within the broad band fil¬ 
ters of Table 1^ to produce mock photometric catalogues. These 
will be used in the next section to assess the reliability of 
AGN/galaxy decomposition using photometric data and quan¬ 
tify the level of accuracy at which physical parameters of the 
underlying galaxy can be derived. Random noise is added to the 
fluxes of each source assuming Gaussian errors with standard 
deviation 10% of the flux in a given waveband. We associate a 
photometric error of 15% with all the flux densities. The broad¬ 
band filters of Table sample a wide range of wavelengths, as 
shown on Figure In order to quantify the level of AGN con¬ 
tribution in each of these filters, we show in Figure the flux 
density ratios of the SEDs that include an AGN component rel¬ 
ative to those that do not (fracACN - 0). A Type 1 AGN with 
a fracACN - 10% will have an AGN emission higher by a fac¬ 
tor of 2 in the UV and MIR rest frame compared to the same 
SED without an AGN component. In UV rest frame, a fraction 
fracAGN - 40% is sufficient to dominate the emission, with an 
AGN contribution four times higher than the young stellar pop¬ 
ulation emission. In MIR rest frame, the AGN emission is domi¬ 
nant for fracAGN >20% with an emission three times higher than 
the stellar emission. The UV and the MIR are thus key domains 
to perform AGN/galaxy decomposition in the case of a Type 
1 AGN SED as also shown in previous work (e.g., [Weedmaii 


[et al.[[2004| [Wu et alT[[2009| l. The AGN emission of the inter¬ 
mediate type is visible in the MIR domains, especially at 4 jum 
rest frame, where it is brighter than the host galaxy by a factor 
of 2 for a fraction of 10%. For the Type 2 AGN, it is clear that 
the emission of the AGN cannot be detected below 2yum rest 
frame. In this model, the torus is optically thick and the emis¬ 
sion from the inner, hotter part of the torus, emitting at shorter 
wavelength is completely absorbed. In FIR rest frame, the emis¬ 
sion of the AGN contributing to 70% of the total F/r will dom¬ 
inate the emission of the host galaxy by a factor of 7. Given the 
ratio between the AGN emission and the host galaxy emission 
(Figure |^, two bands seem to be the key to constrain the Type 
2 AGN emission, the 3-lOjum rest frame and the 30-40/tm rest 


frame as already noticed in previous works (e.g., Faurent et al. 


20001. However, we note that dust emission templates are not 
well constrained in the 30-40 jim range ( [Ciesla et al.|2014 1 and 
thus improving them can help disentangling the AGN contribu¬ 
tion in the FIR. 


3. Recovering the mock galaxy properties 

The SED fitting functions of CIGAFE are applied to the mock 
photometric galaxy catalogue of the previous section to separate 
the stellar emission from the AGN component and investigate 
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Ciesla et al.: Derivation of stellar mass and SFR in AGN host galaxies 
Table 1; Galaxy and AGN parameters adopted to generate mock galaxy SEDs. 


Parameter 

Value 



Description 

Dust attenuation 

E(B - V), 

0.2 




Dust template: 

Draine & Li 

2007 


3.19 Mass fraction of PAH to the total dust mass. 

8.0 Min. intensity of the interstellar radiation field. 

10® Max. intensity of the interstellar radiation field. 

2 Relative contribution between dust heated in photodissociation 

regions, and dust heated by diffuse stellar population. 
AGN emission 


^max ! ^ min 

60 

Ratio between outer and inner radius of the torus. 

Tg.l 

1.0 (for int. type) 

6.0 (for Type 1 & Type 2) 

Optical depth at 9.7 pm. 

P 

-0.5 

Linked to the radial dust distribution in the torus. 

y 

0.0 

Linked to the angular dust distribution in the torus. 


0.001 (for int. type & Type 2) 

89.9 (for Type 1) 

Angle with line of sight. 

e 

100 

Angular opening angle of the torus. 

frac^Qyi 

0., 0.05, 0.1,0.15,0.2, 0.25, 0.3, 

0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.7 

Contribution of the AGN to the total Ljr. 


qpAH (%) 

^min 
^max 

y{%) 


Table 2: Broad-band filter-set used in this paper. 


Telescope/Camera 

Filter Name 

^mean ( 

MOSAIC 

U 

0.358 

HST 

ACS435 

0.431 


ACS606 

0.573 


ACS775 

0.762 


ACS850 

0.9 

Subaru/MOIRCS 

J 

1.2 

CFHT/WIRCam 

Ks 

2.2 

Spitzer 

IRACl 

3.6 


IRAC2 

4.5 


IRAC3 

5.8 


IRAC4 

8 


IRS 16 

16 


MIPSl 

24 


MIPS2 

70 

Herschel 

PACS green 

100 


PACS red 

160 


PSW 

250 


PMW 

350 


PLW 

500 


to assume complex SFHs like those shown in Figure [T] when 
fitting multiwavelength photometric data. It is common practice 
instead, to assume simple functional forms, such as an exponen- 


et al. 2013|), an exponentially increasing SFR (1-exp-ris, e.g.. 

Pforr et al.|2012 Reddy et al.|2012 1, two exponential decreasing 

SFR 

aws with different e-folding times (2T-dec, e.g.. 

Papovich 

et al. 

120011 IBorch et al.||2006l IGawiser et al.| 20071 

^ee et al. 

2009 

', a delayed SFR (e.g., Lee et al.|2010 

2011 

Schaerer et al. 

2013 

', or a lognormal SFH (Gladders et a 

□2^ 

1 . We consider 


in this work the Ir-dec, the 2T-dec, and the delayed models. First 
tests made with CIGALE on the lognormal SFH show that the 
parameters associated are not constrained by broad band pho¬ 


tometry. However, as demonstrated in Figure 2 of Gladders et al. 
( 2013| l, only specific combinations of the parameters of this SFH 
lead to a non null instantaneous SFR at the age of the galaxy. 
Thus, we do not consider here the lognormal SFH as it pro¬ 
vides problematic SFRs. Furthermore, since the 1-exp-ris model 
is recommended to only model the SFH of galaxies at z >2 (e.g., 
Maraston et al.|2()10[ Papovich et al. 201 1[ [Pforr et ah 2012 


Reddy et al.|2012|l, we do not test it in this work. 


The Ir-dec is represented by the following equation: 


how accurately stellar masses and star-formation rates can be 
determined for galaxies that host an AGN. 

To perform the SED htting analysis, CIGALE hrst builds 
models corresponding to a range of input parameters for both 
the stellar and AGN components. The adopted parameters used 
in the htting procedure are presented in Tablej^ The ones related 
to the galaxy host emission templates are selected based on the 
experience gained from galaxy SED modeling at intermediate 
and high redshift using CIGALE (e.g., |Giovannoli et al.||20ir 


Buat et al.|2012[ Burgarella et al.|2013[|Buat et al.|20r4 1 . 

For the building of a galaxy template SED, it is hrst nec¬ 
essary to make some assumptions about the star-formation his¬ 
tory, which will be convolved with the stellar libraries to yield 
galaxy SEDs. The SFH of real galaxies are expected to be highly 
stochastic. It is therefore impractical and probably meaningless 


SFR(t) cc exp{-tlTi) (2) 

where f is the time and ti the e-folding time of the old stellar 
population. The 2T-dec is obtained adding a late burst to the Ir- 
dec SFH, and thus modeled as: 


SFRW^r<-"p „ 

[exp -t/Ti + kx exp -f/ra 


if f < fi - f2 
if f > fi - f2 


(3) 


where T 2 is the e-folding time of the young stellar population, 
and k the amplitude of the second exponential which depends on 
the burst strength parameter fys Finally, the delayed SFH is 
defined as: 

^ fySP is defined as the fraction of stars formed in the second burst 
versus the total stellar mass formed. 
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^restframe (nm) 

lO' 10° lO' 10^ 



^reatframe (/^^) 

lo' 10° lo' 10^ 



^observed 


^reatframe (/^^) 

lo' 10° lo' 10^ 



Fig. 4: The SEDs of Figure normalized to the SED with no 
AGN emission, i.e. fracACN - 0. The three panels correspond 
to one of the Fritz et al. (2006) templates presented in Table 
upper panel is for Type 1 AGN, the middle panel corresponds to 
the intermediate AGN type, and the lower panel is for the Type 
2 AGN. The grey regions indicated different rest frame spec¬ 
tral domains from UV to sub-mm. Normalized SEDs are color- 
coded according to the fracACN parameter. Solid colored lines 
indicate specific SEDs, i.e. where fracACN=^0% (pink), 20% 
(purple), 40% (green), and 70% (red). 


SFR{t) oc fexp(-f/Ti). 


(4) 


We assume a Salpeter IMF and the stellar population models 
of Maraston (2005|l. The metallicity is fixed to the solar one. 


0.02. Dust extinction is modeled assuming the Calzetti et al. 
( |2000| l law with E{B- V), in the range of value shown in Table[^ 
We also assume that old stars have lower extinction compared to 
young stellar populations by a fixed factor, fan = 0.44 ( [Calzetti 
jet al.||2000]l. The UV/optical stellar emission absorbed by dust 


is remitted in the IR assuming the [Dale et al.|(|2014|) tem plates. 
We emphasize that this is different from the Draine & Li| ( 2007} 
libraries used to generate the mock photometric catalogue (see 
Section|^and Table^. This reduces the impact on the results of 
using the same assumptions/templates to generate and to fit the 
photometry of simulated extragalactic sources. 

AGN templates are also included in the fitting procedure. In 
Appendix [A] we demonstrate that in the case of AGN hosts this 
is essential to minimize biases in stellar mass and SFR estimates 
related to contamination of the stellar light by AGN emission. 
Ignoring this effect results in an over-estimation of the stellar 
mass by up to 150% and the SFR by up to 300% depending on 
the spectral type and the strength of the AGN component. The 
parameters used to fit the AGN component are chosen based on 
the results of [Fritz et al.| ( |200^ . However, we decide to fix the 
value of p, y, and 9, that parametrized the density distribution 


of the dust within the torus, with typical values found by Fritz 
et al. ( 2006[l to limit the number of models. Indeed, allowing for 


different values of these parameters would result in degenerated 
model templates. 

Although the AGN parameters used to create the mock SEDs 
are not used in the fitting part of this work, we recognize that 
there could be potential problems with using the same template 
library to generate mock photometric data and then fit them. As 
a test, we have created the mock galaxies with observed AGN 
templates and fit them with the Fritz et aT| ( |2006 i models, the 
results are discussed in Section (3.41 We find that our results and 
conclusions are robust to the set of templates used to model/fit 
the AGN component of the SED. 

The modeled SEDs are integrated into the selected set of fil¬ 
ters, and these modeled flux densities are then compared to the 
ones of the input catalogue of galaxies. For each galaxy in the 
mock photometric catalogue, and for each set of model param¬ 
eters, the is computed. The code then builds the probability 
distribution function (PDF) of the derived parameters of inter¬ 
est (e.g. stellar mass, SFR, fracACN) based on the x^ value of 
the fits. The output value of a parameter is the mean value of 
the PDF, and the associated error is the standard deviation de¬ 
termined from the PDF. We refer the reader to [Noll et al.[P009l l 
for more information on the Bayesian-like analysis performed 
by CIGAFE. 

To avoid any bias that could rise from the fact that we use 
the same tool to create and analyze the mock SEDs, we take the 
following precautions; 


The flux densities of the mock SEDs are perturbed by adding 
a noise randomly taken in a Gaussian distribution with a 
standard deviation of 0.1. 

We use simple analytic SFHs in the SED fitting procedure in 
order to reproduce galform SFHs. 

Although the mock galaxies dust emission is modeled with 
the Draine & Fi ( [2007 [ ) library, we use the Dale et al. ( [2014[ l 
templates in the SED fitting. 

Even though we also use the [Fritz et al.[ ( [200'6l l library to per¬ 
form the fitting, we prevent CIGAFE to use the templates 
used in the building of the mocks catalogues. We did not 
use a different AGN library for our fitting procedure for two 
reasons. The first one is that we are limited by the AGN mod¬ 
els available in CIGAFE. The second one is that it has been 
sh own in [Feltre et al.[ ( [201^ that the smooth torus library 
of Fritz et al. p006[l is highly d egenerated with the clumpy 
torus library of Nenkova et al. (20081. Thus using for in¬ 
stance this very different library will not affect our results. 
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Ciesla et al.: Derivation of stellar mass and SFR in AGN host galaxies 
Table 3: Parameter ranges used in the blind fitting procedure. 


Parameter 

Symbol 

Values 


Star Formation History 

Metallicity 

Z 

0.02 


IMF 


Salpeter 


Double exponentially decreasing 

T of old stellar population models (Gyr) 

Tl 

1,3,5 


Age of old stellar population models (Gyr) 

h 

1,2, 3,4,5 


T of young stellar population models (Gyr) 

D 

10. 


Age of young stellar population models (Gyr) 

h 

0.01,0.03, 0.1, 0.3 


Mass fraction of young stellar population 

fySP 

0.001,0.01,0.1,0.2 


Single exponentially decreasing 

T of stellar population models (Gyr) 

T 

0.5, 1, 3, 5, 10 


Age of stellar population models (Gyr) 

t 

1,2, 3,4,5 


Delayed SFH 

T of Stellar population models (Gyr) 

T 

0.5, 1, 3, 5, 10 


Age (Gyr) 

t 

4, 5, 5.5 



Dust Attenuation 


Colour excess of stellar continuum light for the young population E{B - V), 0.05, 0.1,0.15, 0.2, 0.25, 0.3,0.35, 0.4, 0.5,0.6 

E{B — V), reduction factor between old and young populations /„„ 0.44 


Dust template 


IR power-law slope 

a 

1.5,2, 2.5 


AGN emission 


Ratio of dust torus radii 

^ma.x l^min 

30, 100 


Tgj 

0.3, 3.0, 6.0, 10.0 


P 

-0.5 


7 

0.00 



0.001,50.100, 89.990 

Opening angle of the torus 

e 

100 

Fraction of L/« due to the AGN" 

fracAGN 

-0.2, -0.15, -0.1, -0.05, 0.0, 0.05, 0.1, 0.15, 0.2 
0.25, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8 


Notes. We use low negative values of frac^cN in order to minimize a bias due to PDF analysis, as explained in Section 3.2.1 


We present in Figure [^examples of fits of the mock galaxies 
corresponding to the galform SFH presented in Figure [T] and 
assuming a double exponentially decreasing SFH. In Type 1 and 
intermediate type AGNs, the hts are good but we note some dif- 
hculties in reproducing measurement near the IR peak. This dis¬ 


agreement is attributed to compatibility problems between Dale 
|et al.| ( |2014|l and Draine & Li| ( |2007| l dust emission libraries in 
this range ( Ciesla et al.|2014| l that have no impact on the results 
of this study. For the intermediate type AGN, Figure shows 
that the code uses a model with a strong silicate absorption at 
9.7 yum. However, as shown in Figure (middle panel), the in¬ 
put AGN model used for the intermediate type does not show 
any silicate absorption. Thus, despite the availability of models 
with low values of Tg 7 in the htting procedure, the code does not 
reproduce the absence of silicate absorption of the input SED. 
Although these points have no impact on the result of this study, 
they give a glimpse on possible problems to constrain AGN tem¬ 
plates with broad-band photometry. 


The estimates of the M», SFR, and fracACN as a function of 
the contribution of the AGN are presented in Figure for the 
three different assumptions made on the SFH, as well as the;^^^^^ 
distribution for each case. We discuss these results in the absence 
of AGN in Section 3.1 and the impact of the AGN contribution 
in Section [T 2 I 


Table 4; Fractional differences between the derived M, and SFR 
and the input ones for the mock galaxies with fracACN — 0 (no 
AGN contribution). 


SFH model 

M, 

mean (%) 

o-(%) 

SFR 

mean (%) 

a-(%) 

Ir-dec 

-10.0 

11.5 

-11.4 

7.0 

2T-dec 

-5.2 

10.2 

-6.6 

6.2 

delayed 

-6.5 

17.6 

-9.1 

7.4 


3.1. Determining steiiar masses and star formation rates in 
simuiated gaiaxies without AGN component 

We hrst explore fractional differences between the derived and 
input stellar masses and SFRs for mock galaxies without any 
AGN component included in their SEDs, i.e. the case fracACN = 
0 in Eigure|^ Both parameters depend on the adopted SEH used 
to ht the mock galaxy multi-waveband photometry. Different 
functional forms for the SEH are expected to induce systematic 
variations in both M» and SER determinations ( Bell et al.|2003t 
|Lee et al.|2009||Maraston et al.|2010}|Pforr eFa .|2012| l. Table |4| 

provides the fractional systematic offsets for the different SEH 
assumptions used in our analysis to ht the mock galaxy pho¬ 
tometry. The model that provides on average the best agreement 
between the output and input M* and SERs is the 2T-dec. 

In order to compare the output parameters linked to SEH, we 
show on Figure l^hve galform SFHs as well as the associated 
SFHs obtained from the best ht of the SED htting procedure. 
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Fig. 6 ; Fractional difference between parameters derived from CIGALE and the simulated ones as a function of fracACN^ i e. the 
contribution of the AGN light to the overall IR luminosity. The first column of panels plots the fractional difference in stellar mass 
estimates. The second column is for SFR estimates and the third column corresponds to fracACN- The fourth column of panels shows 
the distribution of Each row corresponds to a set of simulation that assume different input SEDs for the AGN component. Erom 
top to bottom, we present results for Type 1 AGN, intermediate type and Type 2 AGN. Different colors in each panel correspond to 
different functional forms for the SEH. Black is for the 2T-dec model, green corresponds to the Ir-dec model, and purple marks the 
delayed SEH. Dots represent the mean value of the fractional difference at each input fracACN^ the shaded regions show the one cr 
scatter. 


GALFORM SEHs are complex and analytical SEHs cannot repro¬ 
duce the numerous star formation bursts. However, as the inte¬ 
gral of the SEH provide the stellar mass of the galaxy, it is clear 
that, in order to recover the stellar mass, the best Ir-dec and 2 t- 
dec models use a very small, unrealistic age. Indeed, if we fix 
the age of the galaxy, contrary to what is usually done in the 
literature when using these assumptions, by providing a smaller 
age range, between 4 and 5.5 Gyr for our samples as galaxies are 
at z-l, then these two models overestimate by 15 and 19%, 
respectively. Erom the three SEHs models used in this work, the 
delayed SEH better reproduces the global envelop of the gal- 
form SEH. As a result this functional form yields stellar masses 
in reasonable agreement with the input ones and, at the same 
time, a galaxy formation age close to the simulated one. This is 


consistent with the mean SEHs of Illustris presented in Sparre 
|et al.| ( |2014| l. Indeed, the mean and median SEHs of Illustris 
sources have shapes that can be typically modeled with a de¬ 


layed SEH. This result obtained from Semi-Analytical Models of 
galaxy evolution seems to be in agreement with the conclusions 
of Boselli et al. (2001 1 based on observations of local galaxies. 


We thus conclude that a delayed SEH seems to be a more realis¬ 
tic assumption on the SEH of galaxies. 

Several published works have studied the ability of SED fit¬ 
ting techniques to retrieve the stellar mass of galaxies. They were 
based on mock catalogues, on SAMs, or on hydrodynamical 


et al.| 

2009 Lee et al. |2009l Pforr et al. 

120121 IPacifici et al. 

2012 

Mitchell et al.|20131 Buat et al.|2014 

1 . Also using GALFORM 

SEHs, Mitchell et al. (201311 showed that t 

single exponentially 


decreasing SEH provides a good estimation of M» with a small 
offset of -0.03 dex. Using CIGALE, ( |Buat et al.|2014| l found that 
the output stellar mass is systematically lower by 0.07 dex com¬ 
pared to the true one. Despite the different methods used in these 
works, our results are in good agreement as we find an underes- 
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Fig. 7; Comparison between five SFHs from galform and the 
corresponding best-fit output SFHs obtained by CIGALE. The 
green dashed-line is for the Ir-dec model, the purple solid-line 
corresponds to the delayed SFH, and the black dotted-line to the 
2 T-dec model. 


Fig. 5: Examples of hts for mock galaxies associated with one 
of the SEHs shown is Eigure[2 Blue points are the flux densities 
of the mock galaxies, the red lines are the best hts obtained by 
CIGALE. In addition, we show the unattenuated stellar emission 
in black dotted line, the dust emission in black dashed-dotted 
line, and the total emission from the AGN in orange dashed line. 
The associated reduced provided. The single top panel 

shows the best ht in the absence of AGN emission. The two up¬ 
per panels show the best ht for Type 1 with a fraction of 20% 
(left) and a fraction of 60% (right), the middle ones for the inter¬ 
mediate type, and the bottom panels for the Type 2 AGN. 


timation of the stellar mass of ~1% averaged over the three types 
of SEH. 


Small differences are also found in the derivation of the SER 
from the different SEH assumptions, but they globally slightly 
underestimate the SER with offsets between 6.6 and 11.4% 
(Table 1^. These relatively good estimations (<12%) obtained 


with the SER are in perfect agreement with the results of Buat 
|et al.| (|2014[) who found that the SER is robustly estimated. 


with systematic differences lower than 10%, whatever the cho¬ 
sen SEH model as long as one IR data is available. 

Another parameter which is known to be directly linked to 
the estimation of the stellar masses and SER is the amount of at¬ 
tenuation, quantified in this work by E{B-V),. The link between 
the attenuation and the stellar mass is however indirect. Without 
a strong constraint on the dust attenuation, a degeneracy between 
the age of the old stellar population (which is directly linked to 
the stellar mass) and the attenuation appears ( |Pforr et al.|2012t 
Conroy|2013[ Buat et al.|2014| l. We explore this bias by gener¬ 
ating mock galaxy photometry by varying the E{B - V)* input 
value and hxing fracACN-O- The resulting mocks catalogues are 
fit with a 2T-dec SEH model. The systematic offset of the de¬ 
rived M* increases by a factor of two between E(B - V)* = 0.05 
and E(B - V)* = 0.7. The underestimation on the SER varies 
from -15% to -6% for E(B - V)* = 0.05 (Aj/=0.2mag) to 
E(B - V)t = 0.5 (Ay=2mag). Euithermore, if we use, in the 
mock catalogues, a host galaxy SED with a completely obscured 
star formation, then the contribution of the AGN to the total IR 
luminosity is still recovered for high fractions, but the offset be¬ 
tween 15% and 30% is larger by a factor of ~2. Thus the level of 
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stellar light attenuation also plays a role in recovering the stellar 
mass and SFR of galaxies, as already shown in previous studies 


(e.g., Wuyts et al.|2009[ [Mitchell et al.|2013| Buat et al.|2014| l. 


In conclusion, the 2T-dec model provides the best estimates 
of M, and SFR of the simulated galaxies, in the expense of un¬ 
realistic galaxy ages. The delayed SFH recovers the stellar mass 
with a mean offset of ~6.5%, and a mean offset on the SFR of 
~9%, but better reproduces the true SFH of the galaxies. We use 
the results obtained in absence of AGN emission as references 
to analyze the AGN impact on the estimation of M» and SFR. 


3.2. Determining steiiar masses and star formation rates in 
AGN host gaiaxies 

As shown in Figure]^ depending on its intensity, the AGN emis¬ 
sion contaminates large parts of the SED especially the key do¬ 
mains used to retrieve the stellar mass and the SFR of the under¬ 
lying galaxy. In this Section, we discuss the ability of broad band 
SED htting to constrain this contamination, decompose the AGN 
from the host galaxy light, determine stellar masses and SER of 
AGN hosts. The results are presented in Eigurej^for the three 
AGN types considered in this work, where we show the mean 
fractional difference for each parameter as a function of the in¬ 
put power of the AGN (points), as well the 1-cr scatter (shaded 
regions). 


3.2.1. Constraining the AGN contribution 


Eirst we explore the ability of CIGALE to constrain the frac¬ 
tional contribution of AGN emission to the total IR luminosity, 
i.e. the frac^cN parameter. The set of panels in the third column 
of Eigure presents the relative difference between the output 
and input fracACN averaged over the 100 galform SEHs, for the 
different assumptions of SEHs. The AGN contribution is almost 
always overestimated for all three types, except for high fraction 
values ifracAGN >50%) for the intermediate type and Type 2 
AGNs, where it is underestimated. Independent from the input 
AGN SED shape (Type 1, intermediate type, or Type 2), below 
fracACN-^0%, there is a large overestimation of the AGN con¬ 
tribution, that can reach up to ~120%. One explanation could 
be the well-known effect of the use of PDE analysis in retriev¬ 
ing parameters. Eor the lowest value of the parameter, the PDE 
will be truncated, and thus taking its mean value will slightly 
shift the output value of the parameter toward a larger value, 
yielding to its overestimation, as explained in |Noll et al.| ( |2009| l 
and Buat et al. ( 2012| l. However, to prevent this effect, we pro¬ 
vide low negative fracACN as well as input parameters. Without 
any impact on the results of the SED htting and the estimates 
of and SER, it allows a better estimate of very low fracACN 
from the PDE analysis. Thus, the overestimation observed for 
very low fractions shows the difficulty CIGALE has to distin¬ 
guish between a SED without any AGN contamination and an 
SED with a very low contribution of the AGN (~5%). We thus 
conclude that low AGN contributions are very difficult to con¬ 
straint from broad band SED fitting with an overestimation up 
to a factor of 2 for fracACN <10%. At higher fractions, the 
overestimation depends on the type of AGN. A Type 1 contri¬ 
bution to the total L/^ is well recovered with an offset smaller 
than 10% for fractions higher than 10%. In order to better un¬ 
derstand this trend, we show in Eigure|^(Ieft panel) the PDE for 
each fracAGN as seen from the top of it, for the SEDs associated 
to one SEH. We can see from Eigure that the shape and posi¬ 
tion of the PDE do not change when fracAGN increases. Eor the 


intermediate type and Type 2, an overestimation of 30-40% can 
be expected for fracAGN between 10 and 40%. Then the offset 
decreases at higher fractions. These effects are seen on Eigure 
the PDE is large for low values of fracAGN and shows a constant 
overestimation up to high values of fracACN, i-S- 60%. 


3.2.2. Estimating the stellar mass 


The panels in the first column of Eigure |^plot the fractional dif¬ 
ference between the input stellar mass of AGN hosts and that 
inferred from the template fits to the mock photometry. This fig¬ 
ure shows that the stellar masses of AGN hosts can be derived 
with systematic uncertainties smaller than 40% (AlogM* < 
0.15 dex), even in the case of Type 1 AGNs and up to fracAGN - 
0.7. This underlines the importance of AGN/stellar light decom¬ 
position when fitting templates to multi-wavelength photometric 
data of AGN to derive properties of the underlying galaxy. If 
AGN templates were not included in the analysis the inferred 
stellar masses would be biased to high values by a factor as large 
as 2.5 (see Appendix [A|. 


Careful inspection of Eigure]^ also suggests that, in the case 
of Type 1 AGN, the uncertainties on the estimation of M, as a 
function of Type 1 AGN fraction seem to not depend on the as¬ 
sumption made on the SEH up to fracAGN =0.4. The fractional 
difference shows a constant increase up to fracAGN <0.4 and 
then reaches a plateau for higher fractions, except for the 2 t- 
dec model for which the offset decreases. Indeed, as we can see 
from Eigure (upper panel), at an AGN fraction of 20%, the 
AGN emission is higher than the host emission by a factor of 2 
in the end of the NIR domain. As NIR flux densities are known 
to be a proxy for the stellar mass as the emission is dominated by 
the old stellar population (e.g., Gavazzi et al.|1996| , it is sensible 
to think that there is a link between the contribution of the AGN 
in NIR and the variations observed in Eigure To understand 
this trend, we see in Eigure (left panel) the PDE of the stellar 
mass for each fracAGN as seen from the top of it. The PDE is 
very broad and is skewed toward lower values of M* when the 
fracAGN increases. The variation of the fractional difference in 
Type 1 AGNs depends on the contribution of the AGN to the to¬ 
tal Lir. a very weak effect can be attributed to the assumption 
made on the SEH but the Icr scatter of each of them are overlap¬ 
ping showing that this effect is marginal. 


A systematic effect is seen in the intermediate AGN type 
(Eigure first column, middle panel). The offset on the stel¬ 
lar mass slightly increases with the contribution of AGN up 
to fracAGN ~20-30% and then decreases. This threshold of 
fracAGN ~20-30% corresponds to the point where the fracAGN 
starts to be relatively well constrained, even if it is still overes¬ 
timated. We note that the delayed SEH appears relatively less 
affected showing weaker amplitude of variations with fracAGN- 
Given the same pattern observed for three SEHs considered, the 
variation seen here is mostly due to the AGN emission. 


Eor Type 2 sources, the AGN emission has no influence on 
the derivation of M* whatever the SEH chosen. This is likely due 
to the high obscuration of the light emitted from the accreting su- 
petmassive black hole, affecting only slightly the observed rest- 
frame NIR flux of the galaxy. Indeed, it is clear from Eigure ^ 
that the NIR rest frame domain is not contaminated by the AGN 
emission. The M, PDE is relatively narrow and shows the offset 
discussed in Section [XT] (Eigure [9] right panel). 
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Fig. 8; Probability distribution function for the fracACN determination of a particular galform SFH. The vertical axis in all panels 
corresponds to different input fracACN- Each panel corresponds to a different AGN SED added to the stellar light: left is for the 
Type 1 AGN, the middle panel corresponds to the intermediate AGN type, and right panel plots the results for the Type 2 AGN. The 
color indicates the level of probability, blue being the lowest probability and red the highest. The black solid line is the one to one 
relationship. 



Log M. from CIGALE (Mg) 


Log M. from CIGALE (Mg) 


Log M. from CIGALE (Mg) 


Eig. 9: Probability distribution function for the M, determination of a particular galform SEH. The vertical axis in all panels 
corresponds to different input fracACN- Each panel corresponds to a different AGN SED added to the stellar light: left is for the 
Type 1 AGN, the middle panel corresponds to the intermediate AGN type, and right panel plots the results for the Type 2 AGN. 
The color indicates the level of probability, blue being the lowest probability and red the highest. The true log M* for the particular 
simulated galaxy is 10.39 in solar units. 


3.2.3. Estimating the star formation rate 

Eor all three types of AGN considered in this work, the varia¬ 
tion of the SER estimation as a function of AGN strength is the 
same for the Ir-dec, 2T-dec, and delayed SEH. This implies that 
these variations are entirely due to the AGN emission. Eor Type 
1 AGNs, the underestimation of the SER increases with fracACN 
up to 50%. Indeed, as shown on Eigure|^ the UV becomes to¬ 
tally dominated by the AGN emission very quickly. However, 
the EIR is not much affected by the AGN emission which only 
dominates the 30 //m rest frame emission by a factor of 3 for the 
highest fracACN- This allows us to still constrain the attenuation 
and thus the SER, even when there is an AGN contamination and 
estimate the SER with a maximum offset of 30%. The example 
of Eigure (left panel) shows that the PDE is close to the true 
value, but in this case, a secondary peak arises for AGN fractions 
higher than 35%, yielding to the under-estimation. 


In the case of the intermediate AGN type and for input AGN 
fractions of fracACN =10% the host galaxy SER is underesti¬ 
mated by a factor of 20-30%. This is likely related to the over¬ 
estimation of the AGN contribution to the EIR luminosity (see 
Eigure|^third column, middle panel). This results in an underes¬ 
timation of the host galaxy emission in the MIR-EIR domain, 
and thus an underestimation of the SER. However, for input 
fracACN in the interval 25-50%, the AGN contribution to the 
EIR luminosity is better constrained. As a result the inferred 
SER from the SED ht is also in better agreement with the in¬ 
put one. Eor input fracACN ~ 50% the opposite is happening. 
The AGN contribution to the EIR is underestimated and hence 
the host galaxy emission in the rest-frame MIR-EIR part of the 
SED is overestimated. The net effect is an overall overestimation 
of the SER up to 40%. This behavior is also shown in Eigure [T0| 
(middle panel) where the peak of the PDE is slightly shifted to- 
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ward lower values up to fracACN-^0%, then lies on the right 
values up to fracACN—60% where the SFR starts to be overesti¬ 
mated for very high fractions. 

For Type 2 AGNs, we find the same pattern as in the two 
previous cases. The Ir-dec, 2T-dec, and delayed SFH follow 
the same trend. The estimation of the SFR is quite insensitive 
to the AGN contribution with a very weak decreasing up to 
fracAGN—^^%- For higher fractions, the offset on the SFR mea¬ 
surement increases up to an over-estimation up to 20% for the 
2T-dec and the delayed SFH models. A comparison between the 
SED with a AGN contribution of 40% and the SED of normal 
galaxy shows an emission three times higher for the AGN SED 
in the EIR domain, mandatory to have a good estimation of the 
SER ( |Buat et al.|2014| l. The SER PDE corresponding to the 2t- 
dec model (Eigure 10 right panel) shows two peaks, one on the 
right value and one slightly overestimating the SER. Eor high 
AGN fractions, a third peak arises overestimating even more the 
SER and skewing the PDE toward higher values. However, this 
increase is weak for the Ir-dec SEH model yielding to overesti¬ 
mation of 5%-10% at fracACN=T0%- 

3.3. Impact of the photometric coverage 


bands. The absence of EIR and/or submm has no significant im¬ 
pact on the derivation of the stellar mass for any of three AGN 
spectral types. The observed systematic variations with fracACN 
are similar to the case of complete photometric band coverage. 
However, we note a small underestimation of ~10% in M* for 
intermediate type with input fracACN between 30 and 60%. The 
lack of UV photometry has no impact on the determination of M, 
for Type 2 AGN. This is because in this case there is no AGN 
emission to contaminate the UV-to-NIR bands. The absence of 
UV photometry however, does affect the M* estimates for Type 
1 and intermediate type AGN. It leads to an underestimation of 
the stellar mass in Type 1 AGNs, and an overestimation for in¬ 
termediate type AGN for input fracACN = 30 - 60%. This un¬ 
derlines the importance of the UV to recover the stellar mass for 
these types of AGNs. When no data longward of observed 8 pm 
are available, the stellar mass is always moderately underesti¬ 
mated by up to 20-30%, relative to the case of full photometric 
band coverage. This is in agreement with Noll et aL| ( 2009| l who 
found an underestimation of 0.19 dex in the case of local normal 
(i.e. no AGN component) galaxies for the same combination of 
photometric data. 


The results presented in Section ED and Section [T2| have been 
tested in the ideal case where complete photometric coverage, 
from UV to submm rest frame, is available for each source. In 
this section, we study how the lack of photometric data at differ¬ 
ent spectral bands affects the determination of the stellar mass, 
the SER, and the fracACN parameter. As shown in Section ED 
there is no perfect SEH assumption, we thus arbitrarily choose 
the 2T-dec SEH model. This choice does not affect the discus¬ 
sion because we are interested in deviations of the parameters 
inferred from spectral fits using incomplete broad-band coverage 
relative to the ideal case that all spectral bands listed in Table 
are available. We run our SED fitting procedure using exactly the 
same parameters as in the previous sections. However, in succes¬ 
sive runs, the input mock catalogues lack photometry in certain 
broad-band filters. We explore in particular how the results on 
stellar mass, SER, and fracACN change if we exclude in turn, 
the sub-mm rest frame ( Herschel/SPIRE 350 pm and 500 /tm), 
the EIR rest frame (filters with Amem ^ 24 pm in Table |^, the 
mid-IR rest frame (filters with A^^an ^ 8 pm in Table 1^, and the 
UV rest frame (filters with A,„ean < 0.762 pm in Table^. 

Eigure [TT] presents the results for each trial. The impact of 
the different spectral ranges in constraining the AGN contribu¬ 
tion is shown in the third column of Eigure The lack of UV 
rest frame photometry has no impact in the determination of 
fracACN, even for Type 1 AGNs. This parameter is more sen¬ 
sitive to the availability of EIR and submm photometry. Eor all 
input AGN spectral types and for input fracACN ^ 10 - 20%, the 
lack of IR and/or sub-mm data results in a systematic underesti¬ 
mation of inferred fracACN up to typically 20-30% relative to the 
case of complete photometric coverage. Intermediate type and 
Type2 AGN show larger deviations compared to Type 1 AGNs. 
Interestingly, for input fracACN ^ 10 - 20%, the largest offset is 
found in the absence of submm data, showing the importance for 
these long wavelengths in constraining the emission from dust 
heated by young and evolved stellar populations. The exclusion 
of all photometric bands above 8 pm produces systematic offsets 
of up to 50% in the derived fracACN- This is expected because 
CIGALE is based on energy balance and therefore requires long 
wavelength data for optimum performanc e ([Noll et al.|2009| . 

The first column of panels in Eigure |11| shows changes in 
the determination of stellar mass for different set of photometric 


The panels in the second column of Eigure 11 shows that 
the photometric coverage has an impact on the determination of 
SER in AGN dominated galaxies. The strongest effects are ob¬ 
served in absence of data above observed 8 pm. This yields large 
systematic errors of 150-800% in SER, depending on the level 
of input fracAGN- This is in agreement with |Noll et al. ( |20091 l 
who also found an overestimation of the SER by 0.83 dex for 
normal galaxies (i.e. no AGN component) for the same photo¬ 
metric data coverage. Smaller offsets in SER estimates are found 
when the EIR, sub-mm and UV photometric bands are excluded. 
Eor Type 1 and Type 2 AGNs, the absence of any of these spec¬ 
tral domains leads to a systematic overestimation of 5-20%. In 
the case of intermediate type AGNs, the removal of the submm 
leads to variations in the estimation of the SER up to ~20%. The 
situation is worse when no EIR data is available, where the offset 
on the SER increases rapidly with input fracAGN up to 35%. The 
absence of UV data for intermediate type AGNs has a peculiar 
effect as it starts to have an impact from fracAGN-'^0%, where 
the SER is underestimated up to fracAGN-h0%. In intermediate 
types, the AGN has no impact in the UV domain, as shown in 
Eigure]^ and Eigure]^ Removing this range forces CIGALE to 
rely on the EIR to estimate the SER, which is contaminated by 
the AGN emission. It is therefore difficult to constrain the SER of 
objects with AGN SED components similar to the intermediate 
type AGN considered in this work without any UV data. 

Regarding constraints on the stellar mass in the case of SEDs 
without AGN contribution (fracAGN - 0), very small differ¬ 
ences, less than a few percent, are observed when reducing the 
available photometric bands. This is expected since M, is mainly 
constrained from the NIR rest frame emission. The estimation of 
the SER differs by up to ~ 10% when we remove IR and UV data. 
This is to be expected since these spectral domains are important 
to constrain the SER. 


3.4. Impact of the AGN library 

In order to focus our study on the biases that AGN emission has 
on the estimate on M* and SER, we modeled the mock galax¬ 
ies using templates from the Eritz et al. ( 2006| l library, and used 
the same library to perform the SED fitting. Although we are 
cautious not to use the same AGN templates when simulating 
and fitting the mock galaxy photometry, it may be possible that 
the trends we observe are at least partially driven by degenera- 
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Fig. 10; Probability distribution function for the SFR determination of a particular galform SFH. The vertical axis in all panels 
corresponds to different input fracACN- Each panel corresponds to a different AGN SED added to the stellar light: left is for the 
Type 1 AGN, the middle panel corresponds to the intermediate AGN type, and right panel plots the results for the Type 2 AGN. 
The color indicates the level of probability, blue being the lowest probability and red the highest. The true SER for the particular 
simulated galaxy is 3.3 Mg.yr^*. 


cies among the broad-band AGN SEDs of the Fritz et al. 


( |2006l l 


library. We explore this by using a different set of AGN tem¬ 
plates to build the mock galaxy photometric catalogue. We se¬ 
lect three SEDs presented by Lusso et al. ( |2013| l, derived by 
Silva et aT| ( |2004| l, their Seyfert 1, their mildly obscured Seyfert 
2 with log Ah= 21.5 cm“^, and their heavily obscured Seyfert 2 
with log Ah= 24.5 cm“^. These templates bear similarities to the 
the Type 1, Intermediate type, and Type 2 Fritz et al. (|2006|l 


AGN templates defined in Table [T] The approach described in 


Section 2.2 is followed to construct mock galaxy/AGN compos¬ 


ite photometry, with the only difference that the three Lusso et al. 
( |2013| l templates above are used instead of the |Eritz et al.| ( |20()6f 
ones. The SED fitting of the resulting mock galaxy catalogues 


follows the steps described in Section |2i2| (i.e. using Fritz et al. 
|2006| templates). Figure [represents in the same way as Figurejb 
the results on stellar mass, SFR, and fracACN for the 2T-dec 
SFH model. For these parameters the overall systematic trends 
in Figure are similar to those observed in Figure [^ wh ere the 
modeling of the mock galaxies is made with Fritz et al. ( 2006| l 
models. We are therefore confident that our results are insensi¬ 
tive to the adopted AGN template library. 


4. SFR-M* relation 

As mentioned earlier, SED fitting is the most popular method 
to derive the physical properties of very large galaxy samples 
often produced by deep wide area extragalactic surveys. Recent 
studies of such large samples showed that the majority of star¬ 
forming galaxies are known to follow a SER-M* correlation, 
called the main sequence (MS), and galaxies that lie above this 
sequence are experiencing a starburst event ( |Elbaz et al.|201 l| l. 
In order to understand if the results discussed in this work could 
impact the shape of the MS, we show, in Eigure [T^ the mock 
galaxies of our samples in a SER-M* plot. It is known that the 
main sequence predicted by models currently suffers normaliza¬ 


tion problems compared to the observed relation (e.g., |MitcheII 
|et al.||20T4[ [Eurlong et al.||2014| l. At z=l the MS predicted by 
GALFORM is about a factor of 2 lower than the observed MS (see 
Mitchell et al. 2014[ for a complete discussion of this issue). The 
SEHs provided by galform are thus probably not a perfect repre¬ 


sentation of the real z=l star forming galaxies, and it is unclear if 
it could affect our results. Indeed, all the models currently avail¬ 
able suffer from the same issue. However, in this Section, we 
focus on the effect of the AGN emission and the use of SED 
fitting to retrieve the physical parameters of the galaxies on the 
MS. 


As we discuss in Section [3T| and Section |3)2l the choice of 
the SEH assumption has an impact of M, and SER, and thus on 
the SER-M* correlation, as also discussed in||Buat et al.| ( |20T4) l. 
The stellar masses and SER used in Eigure [O] are thos^btained 
from a 2T-dec SEH. The SER-M* relation obtained from the 
Type 1 sample shows a small increase of the dispersion slightly 
shifted toward higher stellar masses, due to the offset on the 
stellar mass obtained for high fractions of AGN discussed in 
Section 3.1 (Eigure [T3] upper left panel). Thus, the increasing 
offset on the stellar mass estimation for Type 1 AGNs has a 
small impact of the SER-M, relation. Eor the intermediate type, 
a shift due to the stellar mass offset is observed, but in addi¬ 
tion there is an increase of the MS dispersion due to the vari¬ 
ation of the offsets with the AGN contribution. Low fractions 
shift the relation toward higher M* whereas the higher fractions 
toward lower masses, thus increasing the dispersion of the re¬ 
lation. However, the Type 2 SER-M, relation does not show a 
significant change as the M* estimation is totally independent 
from fracAGN- Eurthermore, the underestimation observed for 
low AGN fraction in the estimation of the SER seems to not af¬ 
fect the MS of the Type 2 sample. 

The lower panels of Eigure [T^ present the same results but 
for the specific SER (sSER). In contrast with the Type 1 and 
Type 2 samples, the intermediate type sample displays a corre¬ 
lation between the fractional difference between the mock and 
output parameter values and fracACN- The sample without AGN 
is biased toward low sSER and higher stellar mass, and increase 
toward higher sSER and lower stellar mass with fracACN- The 
Type 1 sample only shows an underestimation of the sSER, as 
well as the Type 2 sample where the underestimation is however 
weaker. 

As these results depend on the assumption made on the SEH, 
we provide in Tablej^the linear coefficients of the best fits of the 
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Fig. 11: Evolution of the relative difference between output parameters from CIGALE and the true ones with the input fraction of 
AGN for the three types of AGN considered, normalized to the values obtained with a full photometric coverage. The blue relation 
corresponds to a photometric coverage containing no submm photometry. Green corresponds to the absence of EIR-submm data. 
Orange corresponds to the absence of IR data. Einally, the UV was removed for the magenta relation, distribution are shown for 
every run on the fourth column, we add the distribution corresponding to the total photometric coverage in black for comparison. 
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Eig. 12: Comparison between the output values of M,, SER, and fracACN and the true ones in the case where the AGN emission of 
the mock galaxies is modeled using |Lusso et al.| ( |2013| l templates. The SED fitting is made using the 2T-dec SEH models, and the 
Eritz et al. ( 2006|l models. The results for the Type 1 sample are in blue, for the intermediate type sample are in green, and for the 


Type 2 sample in red. 
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Table 5; Variations on the SFR-M, correlation depending on the 
SFH. 


different conhgurations (i.e. 5 double components and 5 single 
component) for which are derived. 


Coefficients True 

No AGN 

Type 1 

Int. Type 

Type 2 

Ir-dec 

a 

0.68 

0.62 

0.72 

0.47 

0.63 

b 

-6.32 

-5.66 

-6.84 

-4.10 

-5.85 

2r-dec 

a 

0.68 

0.70 

0.74 

0.66 

0.81 

b 

-6.32 

-6.51 

-6.91 

-6.11 

-7.73 

delayed 

a 

0.68 

0.77 

0.56 

0.59 

0.74 

b 

-6.32 

-7.16 

-5.16 

-5.37 

-6.89 


SFR-M* such as: 


log S FR - ax log M* + b. (5) 

We provide these coefficients for the “true” sample, for the sam¬ 
ple containing no AGN, and for the Type 1, Int. Type, and Type 
2 AGN samples, in all of the three SFH cases. In no case were 
we able to recover the true slope of the MS, but the closest value 
is obtained for the 2T-dec SFH (0.70 when the true slope is 0.68) 
which is the SFH assumptions providing the best estimation of 
both and SFR. The largest the offset on the stellar mass esti¬ 
mation, the largest the difference to the true MS slope. It seems 
that the SFR-M, relation is more sensitive to the estimation of 
the stellar mass than of the SFR. The type of AGN has an influ¬ 
ence on the MS slope for all of the three SFHs with a difference 
up to ~32%. 


5. Comparison against iR-oniy SED decomposition 


In this Section, we use our mock catalogues to test another 
method used to disentangle the AGN emission fr om the host 
galaxy in the IR, DECOMPlRp]( |Mullaney et al. 20111. The choice 
of DecompIR is driven by the fact that the code is popular to de¬ 
rive the IR properties of AGN host galaxies and is publicly avail¬ 
able. The aim here is not to compare the results obtained from 
CIGALE and DecompIR as they are two very different meth¬ 
ods. Indeed, results from CIGALE are based on energy balance 
SED htting using multi wavelength data from UV to submm, 
whereas DecompIR is based only on the IR. DecompIR uses a set 
of observed templates for the host and the AGN. The host galaxy 
templates are derived by using the [Brandi et al. ( |2006) l SB sam¬ 
ple, as well as four galaxies taken in the Revised Bright Galaxy 
Sample. These host templates are grouped into hve different av¬ 
erages, spanning a wide range of possible observed SB IR SEDs. 
The AGN template is dehned by effectively subtracting these 


host templates from the IR SEDs of local AGNs (Mullaney et al. 
[ 2 OTTI 1 . 


To determine how successfully DecompIR separates the AGN 
component from the host emission, and thus estimates the SER, 
we perform SED htting with the hve different SB templates and 
an AGN component using MIPS 24 /zm to SPIRE 500 //m data, 
corresponding to 12 to 250 fim rest frame. Eor each mock galaxy, 
we also ht the IR SEDs, using only the hve different SB compo¬ 
nents, i.e. without allowing any AGN contribution. This last test 
is performed to examine the value of adding an AGN compo¬ 
nent to the ht. As a result, each mock galaxy is tested with ten 


The code is publicly available https://sites.google.com/site/ 
decompir/ 


We then use an Akaike Information Criterion (AIC) to select 
the best model, regardless the number of components used to ht 
the data. The AIC is comparable to the more popular Bayesian 
Information Criterion (BIC). Although both criteria are derived 
in the same way, the main difference comes from the prior used 
which is inversely proportional to the number of models for the 
BIC, and a decreasing function of the number of models for the 
AIC. However, some advanced studies comparing both criteria, 
showed that the AIC is more accurate in model selection than 
the BIC (see Burnham & Anderson ( |20()4 1 and [Yang ( 2005 | l for 
detailed comparisons). Eor each mock galaxy, we submit the ten 
different to the AIC. We hnd that for the Type 1, intermediate 
type, and Type 2 AGNs, the AGN components improves the ht 
for 36, 21 and 47% of the galaxies, respectively. Eor all of the 
other galaxies, whatever the type, a SB template is enough. We 
then estimate the AGN fraction, the L/r, and the uncontaminated 
SER, using jKennicuttj ( |T998l l . 

Eigure shows the results, i.e. the variation of the esti¬ 
mation of the AGN contribution, the SER and the L/r with the 
AGN contribution using DecompIR. We present the estimate of 
the L/r as the SER is determined from the simple conversion 
of L/r using Kennicutt ( 1998| l. The AGN contribution is always 
underestimated (except for one point corresponding to fracACN 
~20% in the Type 2 sample). Conhrming what observed with 
the CIGALE results, low fractions of AGN are difhcult to con¬ 
strain. The Type 1 AGN contribution is underestimated by at 
least 25%, whereas intermediate type fracACN is closer to the 
true value with an offset comprised between 15 and 50%. The 
variation of the estimation of the AGN contribution of the Type 
2 sample is however peculiar as it is not monotonic. Like the re¬ 
sults obtained with CIGALE, the SER is slightly underestimated 
but well recovered whatever the fraction of AGN for the Type 1 
and intermediate type sample, and up to 70% for the Type 2 sam¬ 
ple. This is linked to the good recovery of the L/r with a small 
overestimation of a few percent in the absence of AGN, that 
does not evolve with fracACN for the intermediate type sample, 
slightly increases up to ~15% for the Type 1 sample, and follow 
the peculiar behavior observed for the estimation of fracACN for 
the Type 2 sample. However, despite problems in recovering the 
right contribution of the AGN, SER of Type 1 and intermediate 
type AGNs are well recovered. 

The estimate of the SER and L/r of our Type 2 AGN sam¬ 
ple by DecompIR is problematic for fractions larger than 40%. 


One explanation is that Type 2 templates from Eritz et al. (2006 1 


predict IR SED cooler than what was empirically obtained by 
Mullaney et al. (2011 1 . Thus, when DecompIR tries to ht the Type 
2 catalogue created from Eritz et al. ( 2006| ) models, the mod¬ 
els attribute the excess of EIR emission produced by the Type 2 
AGN from Eritz et al. ( 2006| ) to star formation, leading to the off¬ 
set on the L/r that we see in Eigure[^ In order to test this point, 
we run DecompIR on the mock catalogue obtained with the Type 
2 AGN template of jLusso et al.| ( |20T3] l used in Section 3.4 The 
results are presented in purple in Eigure[M The SER and L/r are 
no longer overestimated. DecompIR provides a good estimate of 
the SER with a small underestimation between 10 and 20%, as 
well as for the L/r with an underestimation up to 20% for the 
highest fraction. Interestingly, the estimate on the contribution 
of the AGN to the L/r is well constrained and is similar to what 
obtained with CIGALE, i.e an overestimation of low fractions of 
about 50%. 
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Fig. 13; Top panels; SFR-M* diagram for the Type 1 sample (left), the intermediate type (middle), and the Type 2 sample (right) 
obtained when using the 2T-dec SFH model for the htting. Black stars are the true values of the SFR and M* of the mock galaxies. 
Cyan stars are the galaxies without any AGN component (fracACN-0), purple triangle are galaxies with an AGN contribution lower 
than 30%, and pink circle galaxies with an AGN contribution higher than 30%. Colored-lines are the linear fit corresponding to each 
subsample. Bottom panels; sSFR-M, diagram for the Type 1 sample (left), the intermediate type (middle), and the Type 2 sample 
(right). 



fraCAGN fraCAGN fro-^^AGN 


Fig. 14; Relative difference between the fracACN (left panel) and the SFR (right panel) obtained by DecompIR and the true values. 
Cyan squares, green triangles, and magenta triangles are Type 1, intermediate type, and Type 2 outputs from DecompIR, respectively. 
The purple triangles present the results provided by DecompIR for the Type 2 mock catalogue obtained using the Lusso et ar| ( 20131 l 
templates. The dotted lines indicate the 25% difference. 


6. Conclusions 

We simulate realistic SEDs of Type 1, intermediate type, and 
Type 2 AGNs hosts in order to evaluate the impact of the AGN 
emission on the host SED, and estimate the ability of SED htting 
code to retrieve the physical properties of the galaxy. CIGALE 
is used to model the SEDs from SEHs provided by the SAM 
GALFORM code. Three samples are built, one with the emission 
of a Type lAGN, one with the emission of an intermediate type 
AGN, and one with the emission of a Type 2 AGN, varying their 
contribution through the fracACN parameter. We use the SED 


htting function of the recently updated CIGALE model to derive 
the stellar mass, star formation rate, and contribution of the AGN 
of each mock galaxy, assuming three popular shapes of SEH; Ir- 
dec, 2T-dec, and a delayed SEH. 

In the absence of an AGN contribution, i.e. in normal galax¬ 
ies, the SED htting of our mock samples shows that all the three 
SEHs considered in this work provide good estimations of the 
stellar mass within ~10%. The best estimates are obtained us¬ 
ing the 2T-dec model but at the expense of realistic ages of the 
stellar population, whereas the delayed model provide both. Star 
formation rates are well recovered within 12%. 
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For AGNs, the SED fitting of our mock samples shows that; 

- Stellar masses are overall well recovered with systematics 
up to 40% (0.17 dex) or better, depending on the fracACN 
and spectral shape of the AGN component. This results is 
rather insensitive to photometric bands available as long as 
UV-MIR data are available. is therefore the most robust 
parameter that one can constraint for AGN hosts galaxies via 
broad-band photometric decomposition. 

- The SFR suffers from systematic uncertainties up to 40-50% 
or better as long as FlR/submm data are available. Data-sets 
that are limited to the MIR cannot be used to constrain the 
SFR of AGN hosts. 

- AGN/galaxy decomposition based on broad-band photome¬ 
try can lead to significant overestimation of the AGN fraction 
and hence the inferred AGN luminosity in the case of weak 
AGN. 


We note the need for UV rest frame data to constrain the 
stellar mass and star formation rate of Type 1 and intermediate 
type AGNs. 

The AGN emission has an influence on the slope of the MS 
depending on the SFH assumptions used for the SED fitting and 
the type of AGN . As the AGN contribution can slightly bias the 
estimates of M* and SFR, the variety of AGN types and AGN 
intensities increases the dispersion of the MS. Furthermore, the 
SFR-M, relation is more sensitive to the estimation of the stellar 
mass than of the SFR. 

Finally, we use our mock samples to test a popular method 
used to disentangle the AGN emission from the host emission in 
the IR, DecompIR ( |Mullaney et al.|20TT i, and find that, using the 

“ ' ‘ ( |2006| l AGN templates. 


mock catalogues built from Fritz et al. 


fracAGN is always underestimated but the SFR is well recovered 
for Type 1 and intermediate types of AGN, and overestimated in 
Type 2 AGNs when fracACN >35%. The overestimation of the 
SFR of Type 2 AGNs is due to the FIR prediction of Type 2 AGN 
from [Fritz et al. (|2006| l which is colder than what observed by 
Mullaney et al. l2011|l. When using the Type 2 mock catalogue 


built from jLusso et al.| ( |2013| l, DecompIR recovers well the SFR 
and Ljr with an underestimation up to 20%. 

In this work, we have validated the use of broad-band SED 
fitting methods to derive the stellar mass and SFR of AGN host 
galaxies, as well as the contribution of this AGN to the host 
galaxy SED. This analysis provides the foundation for future 
AGN multi-wavelength studies using CIGALE. Indeed, in the 
near future, the eROSITA X-ray telescope will provide observa¬ 
tions of about three millions of X-ray detected AGNs as well as 
samples of tens of thousands obscured AGNs. These data will 
provide accurate studies of the relationship between accretion 
on the SMBH and the host galaxy properties. Multi-wavelength 
analysis using SED fitting will thus be needed to derive these 
properties. 
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using the 2T-dec SFH model, and not allowing for the use of 
AGN templates. 

Figure A.l presents the fractional difference between the 
output parameters and the true ones as a function of the input 
fracAGN- The estimate of the stellar mass of the Type 2 sam¬ 
ple is not perturbed by the AGN, as we discuss in this work. The 
NIR bands are not affected by the AGN emission, thus the stellar 
mass is well recovered, with or without the use of an AGN com¬ 
ponent. However, Type 1 and intermediate type SEDs are rapidly 
affected by the AGN emission, and M* is overestimated by 50% 
at fracACN=^0%. This overestimation reaches up to 150% for 
fracACN-T0% in the case of Type 1 AGNs. 

The SFR in intermediate type and Type 2 AGNs is overes¬ 
timated by 20% up to /racAGV=30-40%. This overestimation 
increases for higher fractions, especially in Type 2 AGNs where 
it can reach up to 100% The case of Type 1 AGNs is more crit¬ 
ical, the SFR is rapidly strongly overestimated up to 300% at 
/racAGV=70%. 


The xled distributions shows higher value of the ^ 

factor of ~10 compared to the distribution from the htting using 
the AGN component (Figure]^. 

These results show that not taking into account an AGN com¬ 
ponent when performing broad band SED htting of AGN host 
galaxies results in strong biases in the determination of physical 
properties such the stellar mass and the star formation rate. 


Appendix A: The need for an AGN component to 
perform the SED fitting 

In this work, we use a SED htting procedure allowing the possi¬ 
bility of using an AGN component. To understand how omitting 
to use this component would affect the estimate of the stellar 
mass and SER of the host galaxies, we perform a SED htting. 
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Fig. A.l; Evolution of the CIGALE to true relative difference of M, and SER with the fraction of AGN in the case where the AGN 
component is not taken into account in the fitting procedure. Left panel shows the results for the stellar mass, middle panel for the 
SFR, and the right panel present the;^f^^^ distribution. Blue, green, and red data are for Type 1, intermediate type, and Type 2 AGNs, 
respectively. 


20 




















